Tracer Dispersion in Porous Media with Spatial Correlations 
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We analyze the transport properties of a neutral tracer in a carrier fluid flowing through 
percolation-like porous media with spatial correlations. We model convection in the mass trans- 
port process using the velocity field obtained by the numerical solution of the Navier-Stokes and 
continuity equations in the pore space. We show that the resulting statistical properties of the tracer 
can be approximated by a Levy walk model, which is a consequence of the broad distribution of 
velocities plus the existence of spatial correlations in the porous medium. 



The phenomenon of hydrodynamic dispersion — the un- 
steady transport of a neutral tracer in a carrier fluid flow- 
ing through a porous material — has been widely investi- 
gated in the fields of petroleum and chemical engineering 
|Q,D . One can identify different regimes of tracer disper- 
sion according to the Peclet number Pe = v£/D m , which 
is the ratio between the typical time for diffusion £ 2 /D m 
and the typical time for convection ijv. Here v is the ve- 
locity of the carrier fluid, £ a characteristic length scale 
of the porous media, and D m the molecular diffusivity of 
the tracer. 

In the small Peclet number regime, molecular diffusion 
dominates the way in which the tracer samples the flow 
field. In the large Peclet number regime, also called me- 
chanical dispersion, convection effects are significant; the 
tracer velocity is approximately equal to the carrier fluid 
velocity, and molecular diffusion plays little role. The 
tracer samples the disordered medium by following the 
velocity streamlines. In a random walk picture, we may 
think of a tracer particle following the direction of the 
velocity field, and taking steps of length I and duration 
Ijv. 

The classical approach to model dispersion in porous 
media is to consider microscopically disordered and 
macroscopic isotropic and homogeneous porous materi- 
als. Under these conditions, dispersion is said to be 
Gaussian and the phenomenon can be mathematically 
represented in terms of the convection-diffusion equation 
M. This traditional formalism, which is valid for Eu- 
clidean geometries, cannot be adopted to describe the 
global behavior of hydrodynamic dispersion in hetero- 
geneous systems. Specifically, in the case of percola- 
tion porous media, the breakdown of the macroscopic 
convective-diffusion description is a direct consequence 
of the self-similar characteristic of the void space geome- 
try. 

Here we discuss the interesting physics that arises when 
the tracer moves in a flow field with a very broad veloc- 
ity distribution. Consider, e.g., fluid flow in percolation 
clusters near the percolation threshold — a model system 
relevant to a porous medium with stagnant small- velocity 
zones that are linked with large- velocity zones. In this 
case the typical time for convection i/v is without bound 



since the velocity can be arbitrarily small in some fluid 
elements of the void space. Saffman showed |lj that the 
mean square duration of a tracer step is not finite but 
diverges logarithmically unless an upper cut-off is intro- 
duced into the typical time step. This upper cut-off is 
imposed by the mass transport mechanism of molecular 
diffusion. 

Molecular diffusion is expected to affect the tracer mo- 
tion in two ways |Q: 

(i) A quantity of material may cross from one stream- 
line with fluid velocity v to another by lateral diffusion if 
the time step for convection £\\/v is larger than t\, where 
<i = £\/2D m is the characteristic time for molecular 
diffusivity effects to become appreciable || and £\\ and 
£± are the longitudinal and lateral pore lengths, respec- 
tively (with respect to the flow direction) [Q. Thus, if 
£n/v 3> ti, the tracer has enough time to diffuse across 
the pore, and the time step associated with such a move 
is At = ti. When £»/v <C ti, the time duration of a con- 
vective step is smaller than the time required for molec- 
ular diffusion, and the tracer moves with the carrier fluid 
taking a step of duration At = £\\/v. 

(ii) An amount of material may be transported by dif- 
fusion along the pore. The same considerations as in 
point (i) lead to a time step At = £\\/v in which con- 
vection dominates when l^/v <C to = £ 2 /2D m . Here the 
typical length scale is the longitudinal length of the pore 
£u. If £\\/v 3> to, diffusion dominates and the tracer takes 
a time step At = to. 

Here we propose a model of tracer dispersion in a 
porous medium. The porous medium is composed of 
blocks of impermeable material that occupy, with a given 
probability p, a square lattice. We consider a lattice at 
the site percolation threshold, so an incipient spanning 
cluster is formed that connects the two ends of the lattice. 
Previous studies modeled the convective local "bias" for 
the movement of the neutral tracer in the porous media 
assuming Stokes flow Even at macroscopically small 
Reynolds conditions, this assumption might be violated 
in real flow through porous media, specially in the case 
of heterogeneous materials (e.g., percolation- like struc- 
tures) where a broad distribution of pore sizes can lead 
to a broad distribution of local fluxes. 
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FIG. 1. (a) Typical stream-lines of the velocity field in a 
correlated percolation cluster, (b) Velocity magnitudes prob- 
ability distribution averaged over 5 realizations of the per- 
colation clusters, (c) Tracer diffusion in the porous medium 
shown in (a), for Pe = 1.7. We release five walkers and the 
black dots indicate the sites visited at least one time by the 
walkers. 

As a consequence, inertial effects might be locally rel- 
evant. To avoid this problem, we use the steady-state 
velocity field obtained by solving the full set of Navier- 
Stokes in the percolation geometry. Then we study the 
transport properties of a dynamically-neutral tracer mov- 
ing in the flow field. 

We treat the competition between the effects of con- 
vection and diffusion. The velocity field presents a broad 
scale-invariant power-law distribution of magnitude val- 
ues, and we find that there are regions of very small ve- 
locity in which the tracer can be trapped. If convection is 
important, the tracer follows the stream-lines of the fluid. 
When a very small velocity region is reached, molecular 
diffusion effects cannot be neglected, since by diffusion 
the tracer may access the stagnant zones — where it then 



spends a long time. We shall see that due to the existence 
of these stagnant zones, the statistical properties of the 
tracer — e.g., the first-passage time and the root mean 
square displacement — can be understood using a Levy 
walk model for the tracer motion. The existence of Levy 
statistics is also related to the geometrical properties of 
the medium — whether it is correlated or uncorrelated in 
the occupancy variables of the percolation cluster. 

We start by describing the disordered medium and the 
velocity field. Our basic model of a porous medium is a 
percolation model Q modified to introduce correlations 
among the occupancy units || . We assume the existence 
of correlations because we obtain a better mathemati- 
cal representation of transport properties — such as sand- 
stone permeability — by assuming the presence of long- 
range correlations in the permeability fluctuations of the 
porous rock || . The permeability of rocks such as sand- 
stone can fluctuate over short distances, and these fluc- 
tuations significantly affect any fluid flow through the 
rock. Previous models assumed that these fluctuations 
were random and without short-range correlations. How- 
ever, permeability is not the result of a simple random 
process. Geologic processes, such as sand deposition by 
moving water or wind, impose their own kind of correla- 
tions. 

The mathematical approach we apply to describe this 
situation is correlated percolation. In the limit where 
correlations are so small as to be negligible M, a site 
at position r is occupied if the occupancy variable u(r) 
is smaller than the occupation probability < p < 1; 
the variables u(r) are uncorrelated random numbers with 
uniform distribution in the interval [0,1]. To introduce 
long-range power-law correlations among the variables, 
we convolute the uncorrelated variables u{r) with a suit- 
able power law kernel 0, and define a new set of occu- 
pancy variables rj{r) with long-range power-law correla- 
tions that decay as f -7 , where r = \r\ (in the following 
we will set 7 = 0.4). 

We solve the full set of Navier Stokes and continuity 
equations at the percolation threshold of a square lattice 
with 64 x 64 cells and cell edge L = 1 m. Grid ele- 
ment lengths with 1/4 of the solid cell edge, £\\ — £± — 
£ = L/256, have been adopted to discretize the govern- 
ing balance equations within the pore space domain [§. 
Figure |l|a shows a typical velocity field, while Fig. yb 
shows the probability distribution of the velocity mag- 
nitudes averaged over five realizations of the percolation 
clusters. We find that the data are well fit by a broad 
power-law of the type M 
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Next we analyze the transport properties of a neutral 
tracer moving in the fluid. We use a discrete random walk 
model for the tracer motion. Following the Saffman the- 
ory of dispersion in porous media, we define the walker 
motion as a competition between flow-driven convection 
and molecular diffusion. To allow for comparison among 
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different regimes of tracer dispersion, we define a macro- 
scopic Peclet number as Pe = Vi n £/D m , where Vi n = 1 
m/s is the fluid velocity at the inlet boundary of the lat- 
tice. At a given position r in the pore space, we define 
the time scale for convection t c = £/v(r). We choose a 
convective or diffusive move, and a corresponding time 
step At according to: 



t c < td, convection, At = t c 
t c > td, diffusion, At = td- 



(2) 
(3) 



Here td = £ 2 /2D m = Pe£ / (2vi n ) is the characteristic time 
above which diffusion effects become relevant. If the con- 
vection move is accepted, then the tracer moves to the 
nearest-neighbor site in the direction given by the ve- 
locity of the fluid and the clock is updated according to 
t — ► t + t c . If the diffusion move is accepted, then the 
tracer moves to one of the four nearest-neighbor positions 
with equal probability and the clock is updated according 
to t — ► t + td- 

We next discuss the case of large Peclet number, Pe = 
1.7, so the value of td is such that diffusion only occurs in 
regions of small fluid velocity. Typical tracer trajectories 
are shown in Fig. [jjj. We see that the tracer particles per- 
form walks with very long straight trajectories followed 
by periods where they get trapped in small velocity zones. 
These "stagnant zones" in the pore space differ signifi- 
cantly from the dangling ends of the analogous electrical 
problem (i.e., the parts of the infinite cluster connected 
by only one site to the backbone). The tracer enters 
these regions by diffusion, and requires a long time to 
escape. After escaping, the particle performs another al- 
most ballistic trajectory until it penetrates into the next 
small velocity region. The tracer trajectory resembles a 
quasi-one-dimensional channel of "tubes and blobs." The 
"tubes and blobs" picture is the analog for this problem 
of the traditional "links and blobs" picture associated 
with anomalous diffusion in percolation clusters ]s|Jl0[|. 

We analyze the transit time, i.e., the average time re- 
quired for the tracer to traverse a given distance x from 
the inlet line, < x < L, for different Peclet numbers. 
We find (Fig. ||a) that the transit times follow a power 
law 



(4) 



where (3 ~ 1.26 when Pe = 1.7. 

In the "tubes and blobs" picture, we define a tube as 
the set of steps taken by the tracer following a fixed di- 
rection, and we analyze the statistical distribution of the 
tube length s. In stagnant zones where diffusion is dom- 
inant, the tracer is expected to change direction every 
time step, so that s ~ I. In regions in which convec- 
tion dominates, the tracer moves in ballistic trajectories 
limited only by impermeable obstacles. Since long-range 
correlated clusters are very compact (Fig. |l]a), we expect 
s ^> £ and the tube length distribution to be broad. Both 
expectations are corroborated by our calculations. 
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FIG. 2. (a) Transit times for different Peclet numbers, av- 
eraged over five realizations of the percolation clusters, (b) 
Probability distributions of steps for the case of Pe = 1.7 
(power-law distribution) and Pe = 0.1 (Gaussian distribu- 
tion), (c) Transit time exponent as a function of the Peclet 
number. 

Figure ||b shows the probability distribution P(s) for 
two different values of Pe. We find that for sufficiently 
large Pe, the step lengths follow a scale invariant power 
law distribution 

P(s)~s- 2 - 35 [Pc=1.7], (5) 

while for small Pe values, when diffusion is the dominant 
mechanism in the entire pore space, the distribution is 
Gaussian 



p(«) 



= -3( s /so) 2 



[Pc = 0.1], 



(6) 
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with so a characteristic jump length. 

In case (||), the step lengths statistics can be consid- 
ered as a Levy walk, i.e., a random wa lk p rocess in which 
the jump distribution is a power law 

(7) 

A random walker with a distribution (0) travels a typ- 
ical distance r ~ t 2 ~^/ 2 , when 2 < fj, < 3. Thus, the 
transit time for a Levy walker with jump statistics given 

by(|)is0 

(t) ~x 2 '^-»\ (8) 

For fi = 2.35 — the value we find in our simulations for 
Pc = 1.7 — we obtain (t) ~ x 1 - 21 , which agrees with the 
scaling found when we calculate the transit time directly, 
(t) ~ x 1 - 26 from Fig. |2|a, and confirms the validity of the 
Levy walk picture as an accurate description of the tracer 
motion at large Pe. 

The transit time exponent (3 is not universal and de- 
pends on Pe (Fig. |^c). In fact we find that the Levy 
statistics approximates well the value of in the entire 
enhanced diffusion regime 1 < (3 < 2, while in the sub- 
diffusion regime fi > 2, the Levy statistics Eq. (||) ceases 
to be valid. Moreover, we expect two limiting regimes. 
If convection dominates completely (mechanical disper- 
sion), then the tracer should follow the minimum path 
along the spanning percolation cluster. The minimum 
path length £ m i n scales aS ^min 00 ^^min where drnin IS 
the fractal dimension of the minimum path distance be- 
tween two points separated by a linear distance x If 
the tracer moves with a constant velocity, we can iden- 
tify the minimum path distance with the transit time, 
so (3 = d m in- This is the lower limit of the transit time 
exponent, and we confirm this prediction since we obtain 

P ^ ^min 

when Pe is large (Fig. He) Q. 
The other limit at larger diffusivities — the anomalous 
diffusion case |kJ — corresponds to the regime dominated 
completely by diffusion, and the transit time scales as 
(t) ~ x d ™ , where d w is the random walk fractal dimen- 
sion. The value d w depends on the degree of correlation, 
with d w = 2.87 for the uncorrelated percolation limit |5|] 
and d w = 2.41 || for the correlated percolation problem 
we study (7 = 0.4). We see that the limiting cases of our 
calculations agree with these predictions (Fig. ||c). Be- 
tween these two limiting cases, we find that the transit 
time exponent can be approximated by 

/3(Pe) ~ log(Pe). (9) 

We also perform simulations on uncorrelated percola- 
tion clusters. We find a enhanced diffusion regime and a 
sub-diffusion regime as well. However, due to the tortuos- 
ity of the uncorrelated percolation clusters at the thresh- 
old, the distribution of steps is not a scale-free power-law, 



as we find in the case of enhanced diffusion in correlated 
clusters Eq. (jsj) . Thus, we conclude that the Levy statis- 
tics found in the case of dispersion in correlated clusters 
is a by-product of the dynamical properties of the tracer 
moving in a broadly distributed velocity field plus the 
geometrical properties of the particular porous medium 
treated here. The compact features of long-range cor- 
related percolation clusters allows the tracer to perform 
large ballistic steps without encountering obstacles dur- 
ing the random walk process. 

In summary, we find that, at sufficiently large Peclet 
numbers, there is a regime of dispersion for correlated 
porous media in which the trajectory of the tracer parti- 
cle should be better described by a Levy statistics Eq. (JfJ) 
instead of the Gaussian behavior Eq. (||) . Interestingly, 
this fact should be relevant to elucidate the mass and 
momentum transport mechanisms responsible for the dis- 
persion regime called "holdup dispersion" || . Tracer ex- 
periments indicate that this regime of strong dependence 
between dispersion measurements and Peclet number is 
typical of percolation-like porous materials. 
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